Particle entrapment as a feedback effect 
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We consider a suspension of polarizable particles under the action of traveling wave dielectrophore- 
sis (DEP) and focus on particle induced effects. In a situation where the particles are driven by 
the DEP force, but no external forces are exerted on the fluid, the joint motion of the particles can 
induce a steady fluid flow, which leads to particle entrapment. This feedback effect is proven to be 
non-negligible even for small volume concentration of particles. 

PACS numbers: 47.55.Kf, 47.61.-k, 47.15.Fe, 47.54.+r 
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Recent progress and numerous applications in 
medicine, biotechnology and pharmaceutical research 
have witnessed a great interest in understanding of funda- 
mental aspects of particle dynamics in fluid flows at small 
scales [l[ . Particularly, the concept of particle manipula- 
tion is an important problem of micro- and nanofluidics 
Here, the particles can be physical (colloids, liquid 
droplets in microemulsions, small bubbles) or biological 
(cells, bacteria, biopolymers) objects manipulated with 
e.g., electrokinetic H, magnetic forces [4], ultrasound |5[ 
or optical tweezers [6[. Despite a considerable success in 
understanding of the impact of hydrodynamic [3, 7] and 
stochastic [8| forces on the particle dynamics, the prob- 
lem of backward coupling, i.e., the influence of the parti- 
cle dynamics on the fluid, or particle feedback, includes a 
number of open questions. Especially challenging is the 
problem of integral feedback effects that can be induced 
by a collection of jointly moving particles. In this Letter, 
we systematically address this problem. We introduce an 
original theoretical model, apply it to a physically realis- 
tic problem, and predict a generic mechanism of particle 
entrapment that is fundamentally different to the ones 
conventionally known 0, [l(J EH • 

We start with the formulation of a continuum model for 
spherical noninteracting particles of radius a suspended 
in a fluid having viscosity n and density p. To focus 
on the feedback effects we impose a nonuniform external 
force FoF(r) on the particles (Fq its reference value and 
F is the dimensionless field), which, however, does not 
directly influence the carrier phase. Provided that the 
relative size of particles a/L (L is the length scale of the 
flow) and their space- averaged volume fraction $0 are 
small, the dynamics of the two-phase system is described 
by a model allowing for the feedback [121 ]: 
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-Vp + V 2 u + Q s ($)y>F, (1) 
0, v = u + Q s F, (2) 
0, j = ^v-V^ (3) 



where u and v are the fluid and particle velocities, respec- 
tively, p is the pressure, tp is the particle volume fraction, 
and j is the particle flux. The equations have been nondi- 
mensionalized using the scales L, L 2 /D, D/L, r]D/L 2 , 
<3>o (D is the particle diffusivity) for the length, time, 
velocity, pressure, and particle volume fraction, respec- 
tively. Here Q s = 2a 2 LFo/9nD stands for the intensity 
of the external field, Sc = r]/Dp is the Schmidt num- 
ber, and a ratio of two asymptotically small parameters 
($) = 9L 2( 3>o/2a 2 is the feedback parameter, which is 
assumed to be finite. 

To stress the relevance of our approach hereafter we 
provide estimations for a realistic system. As example, 
we stick to the data close to experimental [3] . For a wa- 
ter (rj ~ 10 -2 g/cms, p ~ 1 g/cm 3 ) suspension with 
a ~ 200 nm and <3>o — 3% in a container of the size 
2L ~ 25 pm at temperature 300 K one obtains (<£) ~ 500 
and according to Einstein's formula D 10 -8 cm 2 /s and 
Sc ~ 10 6 . Although in most conventional situations Sc 
is high, it becomes necessary to account for the diffusion 
of particles. The reason is twofold: (i) even small diffu- 
sion gets non-negligible at small scales, e.g., the diffusion 
time is of order L 2 /D ~ 100 s; (ii) it prevents from un- 
bounded physically irrelevant accumulation of particles 
by the external field. 

We emphasize that only particles are able to make fluid 
move, i.e., the fluid flow itself is a perfect indicator for the 
particle feedback, described by the last term in Eq. (pQ). 
Physically, this term comes from the Stokes drag, which 
dominates in the interphase force Io|, EH and is balanced 
by F. Particularly, this results in a distinction of the ve- 
locities of phases, where inertia corrections are negligible. 

Further we focus on an example of a dielectophoretic 
(DEP) force exerted on polarizable particles under ac 
electric field E(r, t) = Re[E(r) exp(iut)]. Hereafter uo is 
the angular frequency, i = Re[z] = z r and Im[z] = 

Zi denote the real and imaginary parts of z. The time- 
averaged force (per unit volume) reads 0: 
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where "*" indicates complex conjugate. The complex 
frequency-dependent function K(uj) = (e p — e m )/(e p + 
2e m ) is a measure of an effective polarizability of the 
particle, known as the Clausius-Mossotti factor. Here, e p 
and e m are the complex permittivities of the particles and 
the fluid medium, respectively. The complex permittivity 
is defined as e = e — ia /uo (e is the permittivity and a 
is the conductivity of the dielectric). The DEP force (|4]) 
comprises two independent contributions: 

F F = Z -e m k r V|E| 2 - Z -e m k % V x (E r x E,) . (5) 

The first term relates to the in-phase component of the 
induced dipole. This force points towards the domains 
of higher field strength for K r > or, conversely, to the 
domains of weaker fields for K r < 0, which is referred to 
as positive-DEP (p-DEP) or negative-DEP (n-DEP), re- 
spectively [14[ . The particles are attracted or repelled by 
the electrode edges. The second term is due to the out-of- 
phase component of the dipole and is essential if there is 
spatially varying phase [3], e.g., for traveling wave DEP; 
it makes the particles move parallel to the electrodes. 

To be able to focus on the feedback effects one has 
to carefully ensure that these are not hindered by other 
possible sources of fluid motion. Because of Joule heat- 
ing, applied electric fields induce temperature gradients 
and therefore create nonuniformities of the conductivity, 
permittivity and density in the fluid, which can lead to 
electroconvection and/or natural convection A typ- 
ical electrothermal force on the fluid Fet <j m UQL~ 3 
is negligible compared to the feedback term, provided 
that Fet/F oc <7 m /7o/<3>o <C 1 (/To is a characteristic 
value of the electric potential). This requirement can 
be safely satisfied for weak fields or weakly conductive 
fluids. Further we assume that the last option is met, 
e.g., for pure water (a m ~ 30 /iS/cm), or deionized wa- 
ter (a m ~ 2 /iS/cm) applied in pjj, and Uo - 0.1 V we 
obtain Fet/F of order 10 -5 or 10 -6 , respectively. As 
we claim below, considerably weaker fields are enough to 
cause the feedback-induced flow (cf. with Uo — 1 — 10V 
m 

HEl EH)- One more advantage of the low conduc- 
tivity is that ac electro-osmosis (another electrohydrody- 
namic effect that can induce a flow caused by electrical 
stresses in the diffuse double layer of charges near the 
electrodes) and natural convection are even weaker than 
electroconvection. 

We now turn to the analysis of a system typical 
for experiments on traveling wave DEP [1, [Hj]. Con- 
sider the two-phase medium filling a rectangular con- 
tainer of sizes L x , L y = 2L, L z and impose a travel- 
ing wave of the potential at the boundaries y = ±L: 
<j) = Uo ex.-p[i(ujt — qx)], where q is the wave number 
and Uo is the amplitude. The complex amplitude 0(r) 
(E = — V(/>) obeys the Laplace equation, V 2 (/> = 0, which 
is readily solved. Assuming that L y <C L x , L y <C L z 
and that the side walls are electrically passive, we obtain 



0(r) = Uo exp(— iqx) (cosh qy)/ coshgL and evaluate 
F(r) = (—Ki coshfo/, sinhfo/, 0) . 



(6) 



Here we define Fo = 3e rn UQq 3 K r /2cosh 2 (b/2) and intro- 
duce a dimensionless parameter Ki = Ki/K r and the 
dimensionless wavenumber b = 2qL. A traveling wave of 
a period 50 jam leads to b « 3 and \Q S \ « 0.5 (estima- 
tions are based on the data as before). Next, we restrict 
our consideration to the case of Ki > as reversal of the 
sign of Ki changes the direction of the induced flow. 

We first point out to a partial case of Ki = 0, which 
corresponds to the limit of perfectly lossless dielectric 
particles. According to (|6]), there is no force allow- 
ing for the particle transport along the plates y = ±1, 
only transversal redistribution occurs. The particles tend 
to migrate either towards or away from these bound- 
aries, which is counterbalanced by diffusion; longitudinal 
nonuniformities are smeared by diffusion. Thus, this sit- 
uation admits a state of mechanical equilibrium described 
by the quiescent fluid u (r) = 0, vanishing particle flux 
jo(r) = 0, and a nonuniform distribution of particles 

,(r) =C *o(j/), *o(2/) =exp(Q s 6- 1 cosh%) , (7) 

where Co = 1/ Jq 1 ^o(y)dy^ as by definition the averaged 
ip is unity. The concentration profile ([7j) describes ac- 
cumulation of particles near the boundaries for Q s > 0, 
or in the center plane for Q s < 0, which corresponds to 
p-DEP or n-DEP, respectively. Analytical and numerical 
treatment of the linearized problem as well as a direct nu- 
merical simulation (DNS) of the nonlinear model ([I])-© 
with (|6]) indicate that the state of mechanical equilibrium 
is stable for any values of the governing parameters. 

What happens in a more general case of lossy particles, 
when Ki ^ and the longitudinal transport is allowed, is 
a simple question to pose, but the remarkably hard one 
to answer. In the limiting case of no feedback, ($) <C 1, 
there is no source for fluid motion and the problem re- 
duces to finding a distribution (/?o( r ) = ^0(^7 2/)? gov- 
erned by Eqs. ([2]), ([3|) with u(r) = 0. In the presence of 
the feedback, the problem becomes highly nontrivial, be- 
cause the system is mechanically nonequilibrium. To get 
an impression of possible scenarios, we have numerically 
integrated Eqs. ([I])-© with dSj) in a two-dimensional rect- 
angular box with the no-slip condition for u and vanish- 
ing normal component for j at the solid walls. A typical 
steady state solution is presented in Fig.[H The flow is of 
a large scale and closed, the particles are involved in the 
vortical motion, which is reminiscent of the particle en- 
trapment und er g ravity by Stommel (for high aspect 



ratio see Ref. 16). However, there are two principal dif- 



ferences. First and most important, the conventional en- 
trapment implies existence of a vortex flow irrespective of 
whether there are any particles or not 0, EH Gil • The 
fluid flow in our system can be induced only by means of 
particles and is not possible otherwise. In contrast to the 



FIG. 1: (Color online). Contourplots of the streamfunction 
(a) and particle concentration (b) in the steady state for 
Sc = 2, Q s = -0.2, Ki = 1, b = 3, ($) = 575. Lighter 
(clockwise rotation for vortices) and darker (counterclockwise 
rotation) colored domains refer to higher and lower values of 
the plotted fields, respectively. Steady states corresponding 
to different Sc look similar: a slight distinction comes from the 
nonlinear term, which is nonvanishing only near the sidewalls. 



previous studies, particle entrapment arises as a generic 
particle feedback effect, which also provides a way to gen- 
erate a flow. The second important distinction is that we 
carefully account for the diffusion effects. This nontriv- 
ial problem has been systematically addressed neither in 
studies cast into a Hamiltonian frame 0, [13, E3] > nor in 
a non-Hamiltionan system [11]. 

A closer inspection of Fig. [1] shows that the flow and 
concentration patterns are one dimensional (ID) every- 
where except for the vicinities of the sidewalls. This al- 
lows for a ID analysis of Eqs. ([I])-© with (|6|) away from 
the sidewalls, valid for systems with a high aspect ra- 
tio. Accordingly, we apply the ansatz tpo(r) = (fo(y), 
uo(r) = (uo(y), 0, 0). As the flow does not influence the 
distribution of particles, (fo(y) is given by ([7j) with the 
same Co. To determine the fluid velocity we account for 
the no-slip conditions at the solid walls uo(±l) = and 
a condition of no mean fluid flux u^dy = 0, implying 
that the flow is closed. As a result, we obtain: 



uo(y) = aV (y), V (y) = V i(y) + V 02 (y) (8) 

with V i(y) = hiX) ~ h(y), V Q2 {y) = -13(1 - y 2 ), where 
a = CoWlQsKil > 0, (3 = 3[J 2 (1) - J 3 (l)]/2 > 0, 
h{y) = J y *o(0 coBhftede, and I l+1 {y) = J Q y I^H (1 = 
1,2); because of symmetry, Vo(y) is an even function. 
Next, we impose the condition of particle entrapment 
J_iJo • e x dy = 0, e x — (1,0,0), which ensures no mean 
particle flux [16[. We arrive at 



(*) = h(l) (c J V *odyj =($) c 



(9) 



representing a formal restriction on ($): for every set of 
governing parameters, only a specific number of particles 
defined by (|9|) can be maintained trapped by the flow. In 
practice, however, e.g., for a closed rectangular box, such 
a restriction is not stringent. For ($) different from ($) , 
the solution away from the sidewalls still corresponds to 
(|7|), as if ($) = (3>) c . Their actual distinction is balanced 
in the vicinities of the sidewalls, where the lack or excess 
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FIG. 2: (Color online). Characteristics of ID state of entrap- 
ment ©-(H} for 6 = 3: velocity profiles (a), maximal absolute 
velocity V m = max y |Vb(y)| and ($) as functions of Q s (b). 



of particles emerges leading to local gradients of concen- 
tration on top of ([7j). This is clearly seen in Fig.QJb). Of 
special attention is the case of ($) considerably smaller 
than (3>) c . Here, although there are not enough particles 
to excite the flow in the whole domain, the entrapment 
still occurs. As before, it is accompanied by the birth 
of a steady vortex flow, but of a smaller longitudinal ex- 
tension. With the decrease of ($) the vortices gradually 
shrink and in the limit of ($) <C 1 no longer exist. 

The characteristics of the discussed state are presented 
in Fig. [21 The particles tend to move along the axis x, 
faster near the boundaries and slower at the center of 
the channel [see Eqs. ([2]) and (SJ)]. Because of viscous 
drag, the fluid is towed by the particles, which has two 
consequences. First, this motion contributes in a positive 
fluid flux, defined by Voi(y). Second, it creates a longitu- 
dinal gradient of pressure that gives rise to an opposite 
Poiseuille flow Vo2(y) with the maximal velocity at the 
center. The velocity profile Vo(y) is a superposition of 
these counterflows, such that the net fluid flux is vanish- 
ing. Noteworthy, profile (|BJ) qualitatively resembles the 
one in the convective flow in a vertical slot induced by 
internal sources of heat [13] • In our case, the role of the 
heat sources is played by the nonuniform distribution of 
particles and the DEP force instead of gravity. Note, the 
state of entrapment exists only for Q c < Q s < 0, see 
Fig.EJb). Beyond this range, Eq. (|9]) prescribes negative 
values of (3>) c , which is physically irrelevant. Particu- 
larly, this manifests that the entrapment can be observed 
only for n-DEP. From the experimental point of view, the 
flow can be easily controlled by tuning the frequency of 
the imposed traveling wave. 

Next question concerns existence of the revealed ef- 
fect in real systems. To answer it, we have linearized 
Eqs. ([I])-® near solution (JJ), (|HJ) and investigated its sta- 
bility with respect to perturbations of the form f(r) = 
f(y) exp(At — ik x x — ik z z). Here A = A r + i\i is the 
complex growth rate, k x and k z are the wave numbers 
along axes x and z, respectively. The analysis has shown 
that the modes with the largest A r correspond to the 
perturbations in the form of rolls, k z = 0. We have 
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FIG. 3: Stability map for b = 3, Sc = 2 (a). Scaling functions 
KxiQs) and K 2 (Q S ) for different Sc (b). 

checked a wide range of Schmidt numbers 1 < Sc < 10 6 
and found no qualitative changes, the stability maps are 
presented in Fig. [3l The regions of stable and unstable 
behavior are separated by two curves k\ c \Q s ) of neu- 
tral stability, on which A r = 0. These lines refer to a pair 
of competing modes of the largest A r and have different 
asymptotes for Sc ^> 1. For the branches with higher 
and lower \Q S \ the scaling laws are K[ c ^ = Kiy/Sc and 
K\ c ^ = K2 Sc, respectively. The dependencies K\(Q 3 ) 
and K 2 (Qs) for different Sc are plotted in Fig. OJb). With 
the growth of Sc, these dependencies converge to master 
curves. Because of the different scaling, the convergence 
of Ki is slower, whereas K 2 gets indistinguishably close 
to its master curve already at Sc = 100. 

The results of the linear stability analysis were con- 
firmed by DNS of the nonlinear model ([I])-® with (|6]). 
Note, as for large Sc the value is high, the instabil- 
ity in this particular situation is hardly reachable exper- 
imentally. However, it can be the case for moderate Sc. 
We have studied the breakdown of the ID state, which 
is found to occur supercritically. The patterns beyond 
the threshold (see Fig. [H for a snapshot) travel along the 
axis x with a speed Xi/k x . For the different branches the 
patterns look similar, but have distinct spatial periods 
and travel in opposite directions. 

In conclusion, we have studied the role of the parti- 
cle feedback in a two-phase system under the action of 
traveling wave DEP. In a situation where the particles 
are driven by the DEP force but no external forces are 
exerted on the fluid, the joint motion of the particles can 
induce a steady fluid flow, which is accompanied by novel 
particle entrapment. In a contrast to the conventional 
mechanism, diffusion of the particles becomes a necessary 
ingredient for the entrapment. This particle feedback ef- 
fect has been proven to be non-negligible even for small 
volume concentration of particles. We note that similar 
phenomena are expected to exist in various physical sys- 
tems. Indeed, the set of Eqs. (PQ)-(J3]) with the force in the 
form (|2j) describe a wide class of problems, e.g., magne- 



tized ferrofluids [l8[ , particles driven by optical tweezers 
@, bubbly fluids under vibrations [19], where the field E 
entering (|4j) is of different nature. 
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FIG. 4: (Color online). Breakdown of ID state beyond the 
stability threshold, Sc = 2, Qs = -0.2, Ki = 7, b = 3, ($) = 
450: streamfunction (a) and particle concentration (b). 

We acknowledge fruitful discussions with A. A. Nepom- 
nyashchy, A. Pikovsky, B. L. Smorodin, V. Steinberg, and 
C. Pooley. S.S. thanks DA AD for support; A.S. was sup- 
ported by the German Science Foundation (DFG, SPP 
1164 "Nano- and microfluidics," project 1021/1-1). 



[1] G. M. Whitesides, Nature (London) 442, 368 (2006). 
[2] T. M. Squires and S. R. Quake, Rev. Mod. Phys. 77, 977 
(2005). 

[3] H. Morgan and N. G. Green, AC Electrokinetics: Col- 
loids and Nanoparticles (Research Studies Press, Bal- 
dock, 2003); A. Ramos, H. Morgan, N. G. Green, and 
A. Castellanos, J. Phys. D 31, 2338 (1998). 

[4] M. A. M. Gijs, Microfluidics & Nanofluidics 1, 22 (2004). 

[5] J. J. Hawkes et al. J. Phys. D 31, 1673 (1998). 

[6] D. G. Grier, Nature (London) 424, 810 (2003). 

[7] S. Kim and S. J. Karrila, Microhydrodynamics: Princi- 
ples and Selected Applications (Butterworth-Heinemann, 
Boston, 1991). 

[8] P. Reimann and P. Hanggi, Appl. Phys. A 75, 169 (2002); 

P. Reimann, Phys. Rep. 361, 57 (2002). 
[9] H. Stommel, J. Mar. Res. 8, 24 (1949). 
[10] M. R. Maxey, Phil. Trans. R. Soc. London A 333, 289 

(1991). 

[11] I. Tuval et al. Phys. Rev. Lett. 95, 236002 (2005). 

[12] Details will be published elsewhere; for approach see 

Refs. [H, El, [H 
[13] O. A. Druzhinin, J. Fluid Mech. 297, 49 (1995). 
[14] H. A. Pohl, Dielectrophoresis (Cambridge University 

Press, Cambridge, England, 1978). 
[15] M. Felten, P. Geggier, M. Jager, and C. Duschl, Phys. 

Fluids 18, 051707 (2006). 
[16] D. V. Lyubimov, A. V. Straube, and T. P. Lyubimova, 

Phys. Fluids 17, 063302 (2005). 
[17] G. Z. Gershuni and E. M. Zhukhovitsky, Convective Sta- 
bility of Incompressible Fluid (Keter, Jerusalem, 1976). 
[18] M. I. Shliomis and B. L. Smorodin, J. Magn. Magn. 

Mater. 252, 197 (2002). 
[19] A. V. Straube, D. V. Lyubimov, and S. V. Shklyaev, 

Phys. Fluids 18, 053303 (2006). 



